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Because of their cosmological origin, gamma-ray burst (GRB) optical afterglows are attenuated 
when they pass intergalactic absorbers in the GRB line-of-sight. Without the knowledge of the 
q number of absorbers and their physical properties, the effect of absorption on the observed mag- 

nitudes can not be determined precisely. Different methods have been applied in order to correct 
for this effect statistically, either using semi-analytical calculations or numerical simulations. We 
follow these works and present the expected magnitude corrections as a function of redshift for a 

t— I set of filters most commonly used in the scientific community. The results are publically available 

1^ on the web (http://igmac.fmf.uni-lj.si). 
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1. Introduction 



Spectra of distant emitters (galaxies, quasars, gamma-ray burst afterglows) reveal a rich set 
of absorption features - a result of light passing through a number of intergalactic gaseous clouds 
(intergalactic medium - IGM). While various absorption features can be found in the observed 
spectra (e.g., Mgll, CIV, OVI, Hell and others), by far the most abundant absorber is neutral hydro- 
gen. Depending on the physical properties of clouds (their redshift, hydrogen column density and 
temperature), the absorption features take various forms (e.g., Lyman forest, Lyman limit systems, 
damped Lyman systems). For a review on IGM absorbers, see [1]. 

Neglecting absorption from all other elements except hydrogen, the light from a distant emitter 
is partially absorbed blueward of the redshifted Lya line. Any photometric observation in these 
wavelengths results in an uncertainty: how much light has been lost due to the IGM attenuation? 

Without simultaneous spectroscopic and photometric observations the photometric correction 
can only be applied statistically, either with the help of semi-analytical tools or simulations. Even 
though the attenuation values obtained in that way are based on average properties of IGM and 
may differ substantially from actual losses of light in some cases (e.g., if a damped Lyman system 
is in the line-of-sight, the attenuation is going to be underestimated), the results are still useful and 
informative. 



2. Methods 



In order to determine the IGM effects on light emitted from a distant source, two methods are 
employed. 

The most widely used method is presented in [2], where measured properties of the absorbers 
(i.e., number evolution with redshift, column density distribution, temperature of the clouds) were 
used to construct a semi-analytical model of the transmission function averaged over many different 
lines-of sight - ^(A). When the transmisson function is known, the attenuation of light observed 
through a filter with known filter transmission function T(X) can be computed as: 

Amag = -1.0861n<2, (2.1) 

where Q is the mean transmission function, averaged over a filter bandpass: 

U ~ $<t>{X)dX ■ K } 

The approach relies on the assumption that the observed flux from a distant emitter at redshift z e m 
is a product of the emitted flux and the mean transmission function. 

The other approach mimics the measurement process. The photometric attenuation is first 
calculated for each line-of-sight separately and then the statistical analysis is applied to the obtained 
results. This can be accomplished with a Monte Carlo simulation and has been applied in many 
works ([3], [4], [5], [6]). Different works mainly differ in the choice of parametrisation describing 
IGM properties. Absorption properties in a particular line-of-sight are determined by a set of 
parameters: 

{zi,N m ,i,bi}i=i,...,N abs , (2.3) 
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where N a b S is the number of absorbing clouds in the line-of-sight, each described by a unique 
set of parameters: redshift z\, neutral hydrogen column density Afai,i and Doppler parameter b\. 
The idea is to simulate a number of lines-of-sight (up to redshift z em ), where iV a b s , Zi, Afai,i and 
b\ are generated according to empirical distributions. Magnitude increment (corresponding to a 
given filter transmission curve) is calculated for each line-of-sight. The resulting distribution of 
increments is then statistically analysed (see below). 

We perform the Monte Carlo simulation by using two different models describing IGM prop- 
erties: one is the same as the one used in [2] (model A), while the other is taken from Table 1 in [6] 
(model B). The latter model is more detailed than the former: (/) it includes two observed breaks in 
the absorbers' number evolution with redshift (e.g., [7] [8]), (ii) it includes a break in the column 
density distribution (i.e., larger fraction of strong absorbers) [6] and (Hi) the Doppler parameter is 
assumed to be distributed according to eq. 3 in [9]. 



3. Results 



We compute the photometric corrections for the redshift interval z em = 0.5 - 6.0 with a reso- 
lution of diem = 0.1 for 10000 lines-of-sight. IGM absorbers are known to have a wide range of 
column density values (from 10 12 — 10 22 cm~ 2 ) [1]. Attenuation for different lines-of-sight there- 
fore differs a lot. The distribution of magnitude increments for an emitter at specific redshift is 
not normal in general (as it turns out, taking the logarithm of the magnitudes does not result in 
a log-normal distribution). Therefore we compute a 50% (± 34.13%) quantile (corresponding to 
median ± 1 a) for each distribution. 

Computed photometric corrections for the case of Bessell U and Bessell R filters are presented 
in Figure 1. The graphs show computed median magnitude correction as a function of redshift 
for: semi-analytical model (blue), simulation results using model A (red) and simulation results 
using model B (green). It is interesting to note that the results obtained by using semi-analytical 
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Figure 1: Magnitude correction for Bessell U and R filters as a function of emitter's redshift: semi-analytical 
model (blue), model A (red; shaded area corresponds to the area between ± 1 a values), model B (green; 
dashed lines correspond to ± 1 a values). Filter transmission curve is given for reference (grey line). The 
redshift and wavelength scales are completely independent! 
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model and model B differ, even though the assumed parametrisation is the same - the difference 
has been also observed in previous works [4] [5]. From the two figures one can notice a large 
uncertainty in the case of Bessell U filter at high redshifts compared to the uncertainty in corrections 
for Bessell R filter. Closer examination reveals that the uncertainty starts rising when most of the 
filter transmission resides below the redshifted Lyman edge (i.e., (1 +z e m) X 91.3 nm), where 
the strength of continuous opacity for high energy photons is highly dependent on the Nm of the 
absorber (e.g., [1] [2] [5]) - hence the large differences among different lines-of-sight. In addition 
to wavelength coverage, amplitude of the uncertainty also depends on the width and shape of the 
filter transmission function. 

4. Summary 

We compute photometric collections for a set of filters most commonly used in the scientific 
community. The results are publically available on the IGM Attenuation Correction web page 1 
and are intended to serve as a quick check on the expected average attenuation of light for distant 
objects at known redshifts. 
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